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Abstract 

We investigated the disentanglement dynamics of two-qubit system in Non-Markovian approach. 
We showed that only the couple strength with the environment near to or less than fine-structure 
constant 1/137, entanglement appear exponential decay for a certain class of two-qubit entangled 
state. While the coupling between qubit and the environment is much larger, system always 
appears the sudden-death of entanglement even in the vacuum environment. 
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I. INTRODUCTION 



A multipartite quantum system, in addition to local quantum coherence that exists within 
each of subsystems, may have nonlocal or distributed quantum coherence that exists among 
several distinct subsystems. This property is entangle, which is superposition of the internal 
states of the systems and cannot be separated into product states of the individual subsys- 
tem. It is recognized as entirely quantum-mechanical effect and have played a crucial role in 
practical application ranging from quantum information^ L 2 L cryptography [2I and quantum 
computation 0, to atomic and molecular spectroscopy 6, 7|. 

Recent, many groups were able to prepare entangled states in a variety of physical systems 
and experimental setups, demonstratin g an impressive ability to manipulate and detect 
them efficiently 
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151 ] In particular, Almeida et al. 15| showed that, 



using an all-optical experimental setup, even when the environment-induced decay of each 
system is asymptotic, quantum entanglement may appear "entanglement sudden death", 
called ESD 18}], that is, entanglement terminates completely after a finite interval, without 
a smoothly diminishing long-time tail. 

As we known, a major obstacle for the controlled entanglement of more and more sub- 
systems remains with the capacity of achieving perfect screening of the system from the 
environment. After some time, the unavoidable residual interaction with the environment 
induces mixing of the system state, and thus the emergence of classical correlations at the 
expense of quantum entanglement. Hence, we face the high relevant task of understanding 
the sources of entanglement decay, what implies the identification of the associated time 
scales. In addition to be enlightened by the experimental discovery of ESD, a lar ge n umber 



.[3,0, 



of theoretical literature have investigated the disentanglement dynamics. [1( % 1171, 118|, |19[ Z. 
Ficek and R. Tanas [16( propose the review with an overview of the mathematical apparatus 
necessary for describing the interaction of atoms with the electromagnetic field. Present the 
master equation technique and describe a more general formalism based on the quantum 
jump approach. F. Mintert et a/.[l3] start with a short recollection of environment models 
adapted for decoherence process in a typical quantum optical context under the assumption 
of complete positivity and Markovian dynamics in the Lindblad form. Yu and Eberly jig , 
showed that the dynamics of the quantum entanglement between two qubits interacting in- 
dependently with either quantum noise or classical noise displays a completely different 
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behavior from the dynamics of the local coherence. Instead of the exponential decay in time 
of the local coherence, quantum entanglement may disappear within a finite time in the 
dynamical evolution. The "entanglement sudden death" has been experimentally demon- 
strated by Almeida. However, it is surprising that few if any fundamental treatment exist 
of decoherence that include the dynamics of disentanglement on better than an Markovian 
approximation or phenomenological. Although, the use of the Markovian approximation is 
justified in a large variety of quantum optical experiments where entanglement has been 
produced, one should notice that Non-Markovian effects are important in the description of 
some condensed-matter system |20l], such as the quantum dot qubit(s) system. Therefore, a 
Non-Markovian effects of the decoherence in any viable realization of qubits is desirable. 

In this paper we examined the disentanglement dynamics of two entangled qubits due to 
spontaneous emission, where the interaction with the environment without rotating-wave ap- 
proximation and the treatment process without Markovian approximation. It is found that 
disentanglement always take only a finite-time to be completed, called "entanglement sud- 
den death" , when the coupling between qubit and environment is strong. While the coupling 
with dissipation environments is weak to fine structure constant 1/137, the disentanglement 
change from exponential decay to entanglement sudden-death with the increasing of the 
portion of the double excitation component in the initial entangled state. We describe the 
sudden-death time of entanglement or the realized lifetime of the given two-qubit entangle- 
ment system through the measurable parameters: coupling constant with the enviornment 
a, energy splitting A and cut-off frequency u c . If we consider entanglements as the central 
resource of most types of quantum information processing, it is the most relevant question 
in entanglement experiment under the environment-induced mixing. 

The paper is organized as follows: In sec. [Til we introduce the Hamiltonian without 
rotating-wave approximation in the two-qubit environment interaction and solve it in terms 
of Non-Markovian treatment. The dependence of the concurrence on the different initial 
condition and the coupling strength to the dissipation environment, are discussed in sec. 
IIHI Finally, the conclusion is given in sec. [TV] 
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II. THE MODEL AND THEORY 



This paper is concerned primarily with two two-level systems, since it is generally believed 
that entanglement of only two microscopic quantum systems (qubits, atoms) is essential 
to implement quantum protocols such as quantum computation. We consider two two- 
level subsystem A, B and assume that each subsystem interacts independently with the 
environment, a well justified assumption wherever the particles composing your system are 
sufficiently separated from each other, and therefore, no collective environment effects must 
be taken into account. In non-rotating wave form such a model may be formulated to the 
following total Hamiltonian (set h = 1): 



H(t) = H qu + H env + H int , (1) 



with 



H qu = -~A A ai-±A B a*, (2) 

H env = ^2 Uka k a k + E U kK b k, (3) 
k k 

Hint = - 9k(4 + a k )a* + - ]T fa(b+ + b k )a* , (4) 

k k 

where the Hamiltonian of the two qubits H qu , the two independently environments H env , 
the interaction H int . Here cr, (i = x,y, z) denotes the usually Pauli spin matrices, Aa (A#) 
describes the energy splitting in the A (B) qubit. (b£), a k and u>k (vie) are t ne creation, 
annihilation operator and energy with wave vector k in the A (B) qubit environment, and 
fk are the qubit-environment coupling strength. Yu and Eberly etc. Q, has employed 
the similar model, but the rotating-wave approximation is valid. Two environments are 
completely defined by the spectral density: 

J ( U} )=^gt8{u-u k ). (5) 

k 

We consider the Ohmic bath J(u) = 2au9(u c — u) in this work, where a is the dimensionless 
coupling constant and 6{x) is the usual step function. 

In order to simplify the non rotating-wave term, we apply a canonical transformation, 
H' = exp(s)H exp(— s) with the generator (2 1|: 

s = E - a ^ + E ~ hk)a *- (6) 

k k k k 



Then decompose the transformed Hamiltonian H' into three parts: 

H = H + H l + H 2 , 



(7) 



where the three parts include the analogous form for A and B qubit, 



(8) 



with 



k k 

Kb = ~V B ^ B af + £ Vkh+hk - £ ^£(2 - £ 

fc k k 

As the same style, = if^ A + H' 1B and = H' 2A + -^2B> where 



(9) 
(10) 



(11) 



with 



H 



2,4 



1 . 

~2 



A 



cosh( V —ik{K ~ h)) ~ V 



sinh(£ -&*))- U £ - &*) 



5% 



^ = exp[- £ ^(Z£)\V B = exp[- £ |l(^) 2 ] 



(12) 



(13) 



(14) 



Here erf = cr A ^fa A : #o * s ^ e Hamiltonian of the noninteracting qubits and environment, H[ 
and H 2 are the interaction Hamiltonian in increasing order of the qubit-environment coupling 
strength g k and f k - Comparing H\ to the term of Hi is replaced by the similar rotating- 
wave approximation term in H[, while the qubit-environment coupling strength g k in Hi is 
replaced by g k r] A A A /(uj k + rj A A A ) in As we seen, g k rj A A A /(uj k + r] A A A ) < g k , that is to 
say, the counter-rotating terms decrease the coupling strength with the environment. 
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Approximately to the order of g\ and /f , we write the total Hamiltonian as H' = H' + H[. 
In the interaction picture, 

V'r{t) = V^A A ^a+afexp [i(w k - rj A A A )t] (15) 



+ V^A A ^a,^exp [-i{u k - r, A A A )t] 
1,1, L j 



K 

+ Y jV B A B ^lbia B exp [i{v k - V B A B )t] 

, ^k 

k 

+ V v B A B ^lb k( r B exp [-i(v k - r, B A B )t] . 
k »k 

We consider in general a system denoted by S interacting with a reservoir or environment 
denoted by R. The combined density operator is denoted by psr. The reduced density 
operator for the system p$ is obtained by taking a trace over the reservoir coordinates, i.e., 
Ps = Tr B (psR)- The equation of motion for p SB is given by 

f t Ps R (t) = -t[V I (t),p SB ]. (16) 

After S transformation, 

j t p's R (t) = -i[Vi(t),p' SR ]. (17) 
This equation can be formally integrated, and we obtain 

Psn(t) = Psr(U) - i J [V;(t),p' SR (t)}dt. (18) 

u 

Here ti is an initial time when the interaction starts, supposing tj = 0. On substituting 

p'snit) into Eq.(18), we find the equation of motion 

t 

j t p SR {t) = -i[v;(t),p SR m - f [v'At), [v;(t),p SR (t)]]dt. m 



We now employ the Born approximation 
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22| in which the interaction between the 



qubit system and the environment is suppose to be weak, and there is no back reaction effect 
of the qubits on the enviornment. In this approximation, the state of the environment does 
not change in time, and we can write the density operator p' SR (t) as p' SR (t) = P^^P^O). 
Under this approximation, Eq.(19) simplifies to 

t 

j t Ps(t)p R (0) = -i\/i(t),p g (0)p a (0)] - / [V;(t), [V;(t ),p' s (t')p R (0)]]dt. (20) 



Substituting Vj(t) into Eq.(20) and assuming that the environment modes is in thermaliza- 
tion, the Tr R are given by: 



Tr R [btb k p R ] = Tr R [b k p R b+} = n k , 
Tr R [b k b+p R ] = Tr R [bt PR b k } = n k + 1. 



(21) 
(22) 



Then, 



d i , . 
dt P ^ 



(23) 



$> A A A ^) 2 (n^ + l)[p5(i>^- " ^-Ps(*'K] exp i(u k - V A A A )(t - t) 



dt 



u k 



dt 



f J2(v A &A 9 ^) 2 n£[ps(t)^4 - ay s (t)o A ] exp \-i(u) k - V A A A )(t - t)} dt 



- I ) 2 ("* + !) - *-Ps(t'K] exp -i(u; fc - ^)(t - f' 



v k 



J2(v B *B^) 2 (n» + l)^(f')*?^ " °-Ps(t)^} exp - „ fl A fl )(f - t) 



dt 



v k 



Y,{v B A B f ^fn B k [p n s {t)a B a B + - a B p'(t')a B ] exp -i{v k - rj B A B )(t - t 



dt 



v k 



dt 



;>> BA *— )Vf + 1) W^(t') - <7^(*'K] exp - rj B A R ){t - t) 

, V k L 



In this equation, the n k and n& + 1 term on the right hand side describe, respectively, decay 
and excitation process, with rate which depend on the temperature, here parameterized by 
n k , the average thermal excitation of the reservoir. In this work, we study the limit of zero 
temperature, n k = 0, that is to say only the spontaneous decay term survives leading to 
purely dissipative process. 

The matrix equation is solved in the representation spanned by the standard two-qubit 
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product states basis |1) = | ft)? |2) = | TD? |3) = | IT)? |4) = | II)- After Laplace transfor- 



mation and convo 
obtained as follow 



ution theorem, the master equation of the system of two qubits can be 



23| : 



(r) A A 



Pp's(P) ~ Ps(P) = ~ E P V ' i(ll A ZL M W4<* - °-Ps(P)<\ (24) 

('/"A/f 



_ y W Uh AK , [a A a A p'(P) - a A p'{P)a A ] 
^ P + i{u k - r] A A A ) 1 + ~ P5V ; ~ HSK ' +J 

^ P + i{v k - T) H A B ) 

( V A A A ^-) 2 

Denote the summation of the environment degree of freedom A . = > , 75—7 w \ . > . A_ = 

( J? A AA £ft|i)2 (/A B ^) 2 

E fc p^k-^a,) . 5 + = E fc P +8 n-^A B) aild B - = z2k p- i{ v k -$A B y The deca y rate 18 

dependent on the process, as seen from A + , A_, B + and £L, instead of constant for all 
process in Markovian approximation. We shall therefore focus on the precise time scales of 
every decay process. 

According to the Kronecker product property and technique to Lyapunov matrix equation 
in matrix theory, expand matrix into vector along row of the matrix from two sides of master 
equation, 

|P/16X16 + [A_/ 4X 4 ® ((X+af <g> I 2 x2) T + S_/ 4 X4 ® (hx2 ® <^-) T ] (25) 

-L4_(a A g> I 2x2 ) <g> (a A g> J 2 x2) T + 5_(/ 2 x2 ® ^ B ) ® (^2x2 ® ^+) T ] 
-[A+(<7* <8) J 2X2 ) <g> (a A <g> J 2x2 ) T + fi+(/ 2x2 ® ^ B ) ® (I 2 x2 ® ^+) T ] 
+ L4 + (o^o-f ® 7 2x2 ) ® J 4X 4 + S_(/ 2x2 ® a^<x B ) ® /4x4]|vec[^(P)] = Vec[// S 

The 4x4 matrix equation transformed into 16 x 16 matrix equation with the form 



U(P) 16xl6 Vec[p' s (P)) = Vec[p' s (0)} (26) 



where Vec[p' s (P)] is the vector of row expanding of matrix p' s (P). The solution formally is 



Vec[p s (P)] = U(P)^ xl6 Vec[p' s (0)}. (27) 



Inverse Laplace transformation to time parameter space, 



J^-Vec[p' 5 (P)] = JZ-'UiP^Vecip'siO)}. 



(28) 



i.e. 



Vec[p' s (t)} = Jf-'UiP^Veciip'siO)}}. (29) 

^^ 1 f/(P) 16 1 xl6 can be obtained (see Appendix). 

Compared with Markovian approximation, decoherence rates 7(0;) in our results becomes 
frequency dependent. Due to entanglement and environment interaction together, the decay 
rate for variety process are different, some increase slower, some increase faster, i.e. the two 
bathes has indirect interaction through the two entangled qubits. That is more general and 
physical. 

Therefore, the reduced density matrix p' s (t) in the Schrodinger picture is obtained p' s (t) = 
exp(— iH' t)p" s {t) exp(iH' Q t), the matrix form is 



Ps(t) 



Ps(t) 



exp(i^At) 





exp(— i 







exp(— i- 1 — ^ 



/ exp(?2^t) 



expH^t) 



(30) 







exp( 



■t) 







exp(i^t) / \ exp(i^t) 

Transform p' s (t) into ps(t) through ps(t) = Tr R [exp(—S)p' s (t)p R (0) exp(S)], denoting X A 
E k fttfto ~ a ^ X * = E k t k ^(H - hi so 

ps{t) = Tr R [(coshX A - smhX A a A ) <g> (coshX B - sinhX B crf )p s {t)p R 



(31) 



(coshX^ + smh.X A (j A ) <g> (cosh Ab + sinhX^uf ) 
1 + r] A 1 + rf , l + r] A l -r] B 
— — Ps(t) + — T 

1 - T] A 1 + T] B 



/2x2®af)p 5 (t)(/ 2x2 ®0 



2 2 

1 - 7] A 1 - 7] B 



(a A ® I 2x2 )p s (t)(a : 
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a';®a B ) Ps (t)(a A ®a B ). 



2 2 

Until now, we obtain the reduced density matrix in all kinds of initial state. 

Although a general solution to this problem, for arbitrary system dynamics and deco- 
herence mechanisms is still out of reach, out technical machinery, developed in the previous 
section allows to treat arguably all situations encountered in typical state of the art experi- 
ments, as in the quantum optics and condensed matter. 



III. THE RESULT AND DISCUSSION 



We assume that at t=0, the two qubits and environment are described by the product 
state exp(— S)ps(0)pr(0) exp(S') = ^ab ® |0)a|0)_b, where ^ ab is the entangled initial state 
of the two qubits and |0)yi|0) ^ is the vacuum state of two environments. Let us assume 
that the initial density matrix is only practically coherence of a familiar type (one of the 
atoms is excited, but it is not certain which one). This is easily expressed in the following 



form 



19 



24| 





fd 








o\ 


1 





c 


z 





3 





z* 


b 







V° 








a) 



(32) 



where the factor 1/3 is for notational convenience. In order to compare with previous results, 
consider an important class of mixed state with single parameter a satisfying initially a > 0, 
d = 1 — a, and b = c = z = 1. We will use Wootter's concurrence to quantify the degree of 



entanglement [25l. l26l| . Let p be density matrix of the pair of qubits expressed in the standard 
basis. The concurrence may be calculated explicitly from the density matrix p for qubits A 
and B: C = max(0, \/X\. — s/Xz — \/A~3 — VXi), where the quantities Aj are the eigenvalues of 
the matrix M: M = p{ay ® <jy)p*((jy ® a^), arranged in decreasing order. Here p* denotes 
the complex conjugation of p in the standard basis. It can be shown that the concurrence 
varies from for a disentangled state to C — 1 for a maximally entangled state. 

Firstly consider very weak qubit-environment interaction, oia = olb = 0.01, which is 
larger a bit than the fine-structure constant 1/137. Here and in the following, energies A^, 
and Ab are expressed in units of u c , times in units of uj~ x . We assume A^ = 0.2, A# = 0.4. 
In Fig.l, the time evolution of the concurrence for various values of the parameter a is 
shown. The figure shows that for all a values almost between 0.3 and 1, concurrence decays 



is completed in a finite-time, which is the effect of "entanglement sudden death" [15|, [18], 
but for smaller a's the time for completed decay is infinite, which is consistent with Ref.15 
and 19. The result indicated that in the weak dissipation environment, such as the all- 
optical setup in Ref.15, the Markovian approximation and rotating- wave approximation are 
available. When the coupling constant to the enviornment is near to or less than fine- 
structure constant, we see that the quantum dissipation of the vacuum environment is 



9 



not sufficient to completely destroy the entanglement in a finite time in some situations. 
The sudden death of entanglement results from the decays of the mixed double excitation 
state component. With increasing of the mixed double excitation state component, a value, 
concurrence change from exponential decay to sudden death. The entanglement has another 
unusual relaxation property: different entangled states, corresponding to different values of 
a, with the same initial degrees of entanglement may evolve with different route, some 
showing entanglement sudden death, some not, some decay faster, some slower. That is to 
say, we can prepare certain initial state to prolong entanglement time. 

Next, consider large qubit-environment interaction, a a = olb = 0.05, the other parame- 
ters and initial entangled state are same with Fig.l. The time evolution of the concurrence 
through the entire range of different a values is plotted in Fig.2. As we shown, concur- 
rence actually goes abruptly to zero in a finite time and remains zero thereafter. That 
is to say the entanglement sudden-death always happens. In the first example above, we 
have shown that the entanglement can last for infinite period in the vacuum reservoir for 
some initial entangle state. However, in Fig.2 the sudden death of entanglement always 
happens no matter which entangled state the qubit are initially in. That is also shown 
that the disentanglement dynamics varies with the coupling strength with the enviornment 
or the rotating-wave approximation and Markovian approximation is unavailable, when the 
coupling to the environment is much larger than the fine-structure constant. Fig. 3 shows 
concurrence for a a = olb = 0.1, under the same initial condition. It is observed that in 
the same initial state, the death time decreases as the increasing of the strength of qubit- 
environment interaction. 



IV. CONCLUSION 



In this paper, we consider two two-level qubits that are spatially separated from each other 
and independently coupled to local vacuum environments. We investigated the dynamics 
evolution of entanglement between the qubits. We show that, for a certain class of two-qubit 
entangled state, the entanglement measured by concurrence can change from exponential 
decay to sudden death with increasing of the mixed double excitation state component in 
the case of weak coupling with environment. Increasing coupling strength, the entanglement 
sudden-death always happens no matter which entangled state the qubit are initially in. The 
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exponential decay of entanglement is a very special result to the weak dissipation vacuum 
reservoir. The entangle sudden death time in our result is obtained from the physical 
parameter: coupling constant a, energy splitting A^, A B and cut-off frequency lj c . Finally, 
we hope that this work will stimulate more experimental and theoretical works in quantum 
information and computation for quantum optical control. 
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Appendix 

In this Appendix, we give details of how to inverse Laplace transformation to time 
parameter space. ^(-P);T6xi6 * s composed by the matrix element: p ^ A , p ^ A , 
P+A-+A+' P+A-+B+' p+a_+a + +-B+' p+A-+a\+b++B- etc - Tnen ^ lu ( P )iixi6 is inverse ev- 

p 1. SOxvv ;)J _ i 



ery matrix element. As we know, j£f ^ = 1. Solve Jzf 1 p 1 . etc. through the following 



method. 

cr+ioo 

1 1 f exp(Pt) 



J*?- 1 — ^— = — / ^ ° J A dP (Al) 

P + A_ 2m J („a Aa £^. )2 V ; 



(J — 200 



Then Changing P to ilj + + .[27] 



cr+ioo +oo 

I f eMPt)^ dp= l[ »PN^) ^ (A2) 



2m J {v a Aa ^ )2 2m J ( „a Aa £^ )2 

r ~T~ l^k P~i(u; k - V A A A ) W l^k (u J + V A A A )-u Jk -i0+ 



Denote R(lj) and j(cj) as the real and imaginary parts of ^2 k (r] A AA^j^) 2 / (lj — uj k — i0 + ) 



R{u) = Y,P 



( v a a. 



9ktk\2 



LJ - LJ k 



(r)Afp I dJ 

10 



J{J) 



-2a 



lj c 



lj 



Ml! 

lj + r}A \ lj c + i]A lj + r]A 



In 



(lj-lj')(lj' +T]A) 2 

Lj \ (lJ c + TjA) 



i]A(lj c - lj) 



(A3) 
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and 

= 2a-KU 1 (r?A) ' ■ (A4) 

Where p stands for Cauchy principal value. 

+00 

jy-i 1 = J_ / expM + 0+) ^. fA5l 

P + A_"2m]uj-R(uj + T] A A A ) + i-fiu + r) A A A ) { ' 

—00 

= exp[iu; i^ - 7(^01 + ^ A A A )t] 

where u>o\ is the solution of equation uj — R{lo + T] A A A ) = and is the Lamb shift due to 
the local interaction of the qubit with the enviornment. 
In the same way, 

P \ A = ex P[-^oit - 7(% + r) A A A )t\. (A6) 



It is clear that ££ p . conjugate with Jzf 



p+A- ^ "j "a""- — p+A+- 



-oo 



^ r~ , = — / ?^ + 0¥) , „ A > (AT) 

P + A_ + A+ 2m J cu + ij(r] A A A )+i-f(ri A A A ) 



-oo 



= exp[-2 7 (r / A A A )t], (A8) 

The decay for Jzf -1 p^ 1 ,^ accelerated (by a factor of almost two) as compared to 
-l i 



Jz? p+x~ ' un der the influence of zero temperature environment. [1 

i (A9) 

+00 

1 /" exp(iu;t + 0+) ^ 

" 2m J lu-R(uj + r] A A A ) - R(lu + rj B A B ) + i^(u + i] A A A ) + z 7 (cj + rj B A B ) U 
—00 

= exp[iu s 12 t - + f] A A A )t - + T] B A B )t], 

where u{ 2 is the solution of uj — R(uj + r] A A A ) — R(u + t] b Ab) = and is the Lamb shift 
due to the two environments indirect interaction, which is a nonlocal effect. Jzf -1 5 , . 1 , 5 

7 -r +-D4- 
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conjugates with Jzf 1 P+A l +B . 

r~ — (A10) 

+00 

1 / expficut + 0+) 

-doj 



2ni J uj- R(r) A A A -uj)- R(uj + r) B A B )+i>y(ri A A A -Lj)+i>y(uj + ri B A B ) 

—00 

= exp[iw? 2 f - 7 (^ A A A - u a 12 )t - 7 K 2 + v B ^B)t], 
where uj® 2 is the solution of uj — R(t] A Aa — uj) — R(uj + r] B Ab) = and is also the Lamb shift 



due to the two environment indirect interaction. In the same way, Jzf 1 P+A +B+ conjugates 



with if 



-1 1 



P+A++B- ■ 



; — 7T- (All) 

+00 

— / exp{iut + + )/[u - R(t] A A a -uj) - R{t] A A a + uj) - R{uj + r] B A B ) 
2iri J 



-00 

+i 1 {r t A A A -u) + h(r] A A A + uj) + i-f(uj + 7] B A B )}du (A12) 
= exp[iw 3 it - j(r] A A A - uj 31 )t - -f(r] A A A + uj 31 )t - 7(^31 + ?7 B A B )i], 

where 0*31 is the solution of w - R(r] A A A — uo) — R(r/ A A A + uj) — R(uo + r] B A B ) = 
and is another Lamb shift due to nonlocal interaction. P+A+ + A +B+ conjugates with 
if 



-1 



P+A++A-+B 

9-1 



if — — (A13) 

P + A_ + B + + v ; 

+00 

— / exp(icjt + + )/[cj-i?(cu + ^A j4 ) - R(t] B A b -uj) -R(t] B A b + uj) 
2m J 



-co 

+i-f(u + r] A A A ) + i-f(r] B A B - uj) + ii{r] B A B + u)]du 
= exp[io; 3 2t - 7(^32 + ri A A A )t - ^(i] B A B - uj 32 )t - <y(i] B A B + uj 32 )t], 

where u 32 is the solution of uj — R(uo + r] A A A ) — R(r/ B A B — uj) — R(r/ B A B + uj) = and 
is the Lamb shift due to the two environment indirect interaction, too. J?" 1 p+A + B++B 
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conjugates with 1 p+a+ +b + +B- - 



(A14) 

P + A + + A_ + B + + B_ y ! 



= -f 

2m J 



expUut + + ) 

duj 



oo + i^y(r] A A A — + i^y(r] A A A + u) + ij(u + r] B A B ) + i^(i] B A B - u) 

oo 



exp[-2 7 (77 A A A )t - 2 1 { V B A B )t}. 
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Figure Caption 

Fig. 1. The entanglement decay via spontaneous emission of two two- level qubits starting 
from the initially entangled state (1 - a)/3| |T) < T T I + a/3| II) (III + 1/3(| IT) + | TI 
))((TII + (ITI) with a between zero and 1. the coupling constant of the environment 
and qubit a a = olb = 0.01. Here and in the following figures energies A^ and Ab are 
expressed in units of u c , times in units of We assume A^ = 0.2, A B = 0.4. 

Fig. 2. The entanglement decay via spontaneous emission of two two-level qubits the 
coupling constant of the environment and qubit ctA = olb = 0.05. the other parameter 
the same as Fig. 1. 

Fig. 3. The entanglement decay via spontaneous emission of two two-level qubits. the 
coupling constant of the environment and qubit a a — olb — 0.1. the other parameter 
the same as Fig. I. 
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Fig.2 




Fig.3 




